Proposal for the detection and braiding of Majorana fermions in a quantum spin Hall 

insulator 
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We show how a quantum dot with a ballistic single-channel point contact to a superconductor 
can be created by means of a gate electrode at the edge of a quantum spin Hall insulator (such as 
an InAs/GaSb quantum well). A weak perpendicular magnetic field traps a Majorana zero-mode, 
so that it can be observed in the gate- voltage-averaged differential conductance {dl/dV} as a 4e 2 jh 
zero-bias peak above a (§7r 2 — 4)e 2 //i background. The one-dimensional edge does not permit the 
braiding of pairs of Majorana fermions, but this obstacle can be overcome by coupling opposite 
edges at a constriction, allowing for a demonstration of non-Abelian statistics. 



Topological insulators in proximity to a superconduc- 
tor have been predicted [T to support Majorana zero- 
modes: midgap states with identical creation and an- 
nihilation operators and non-Abelian braiding statistics 
[2 [3] , that are presently under intense scrutiny [4 . The 
conducting edge of a quantum spin Hall (QSH) insulator 
seems like an ideal system to search for these elusive par- 
ticles in a transport experiment 0E]: Only a single mode 
propagates in each direction along the edge, unaffected 
by disorder since backscattering of these helical modes is 
forbidden by time- reversal symmetry [7] . The QSH edge 
is thus immune for the multi-mode and disorder effects 
that complicate the Majorana-fermion interpretation of 
transport experiments in semiconductor nanowires [5] [S] • 

Andreev reflection at a superconducting interface has 
been reported in an InAs/GaSb quantum well 10) , which 
is a QSH insulator because of a band inversion and the 
appearance of edge states connecting conduction and va- 
lence bands |llj . Similar experiments can be tried in 
HgTe/CdTe quantum wells, where the QSH effect was 
first discovered |T2l [13]. We expect a Majorana fermion 
to be present in these systems, delocalized along the edge 
connecting a normal and superconducting contact, but 
without a distinctive resonance in the electrical conduc- 
tance. Andreev reflection of a helical edge mode doubles 
the current at all energies inside the band gap, so each 
edge contributes 2e 2 /h to the differential conductance ir- 
respective of any midgap states. 

Here we present a method to restore the sensitivity of 
the conductance to the zero-mode resonance, by trap- 
ping the Majorana fermion near the superconducting in- 
terface. Only a minor modification of the existing ex- 
perimental setup [TU] is needed, essentially only a gate 
electrode at one of the edges, to locally push the conduc- 
tion band through the Fermi level. (See Fig. [TJ) The area 
under the gate then forms a two-dimensional metallic re- 
gion, connected to the superconductor by the helical edge 
mode. Backscattering at this Andreev quantum dot in 
a weak magnetic field (one flux quantum or less through 
the dot) provides for an electrostatically tunable confine- 
ment of Majorana fermions. We discuss the detection of 
Majoranas as a short-term application, and braiding as 
a longer term perspective. 
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FIG. 1: Left panel: Andreev quantum dot, created by a gate 
electrode at the edge of a quantum spin Hall (QSH) insulator 
in a perpendicular magnetic field B. A current / is passed 
between metallic and superconducting contacts, and the dif- 
ferential conductance G — dl/dV is measured as a function 
of the bias voltage V for different gate voltages. Right panel: 
Band structure of an InAs/GaSb quantum well, for the pa- 
rameters used in the computer simulations. The helical edge 
states appear inside the gap, connecting conduction and va- 
lence bands. Below the gate, the conduction band is pushed 
through the Fermi level at Ef, to create a metallic puddle. 
Inside the superconducting contact, a gap 2A opens around 
the Fermi level. 



There exists a variety of phase coherent backscatter- 
ing mechanisms for helical edge modes [14—20], based on 
different methods of time-reversal symmetry breaking to 
open a minigap in the edge state spectrum. A locally 
opened minigap forms a tunnel barrier for the edge modes 
and two tunnel barriers in series form a quantum dot at 
the QSH edge [19]. For a robust Majorana resonance it 
is advantageous to have a ballistic coupling rather than 
a tunnel coupling to the superconductor, so we form a 
quantum dot by placing two ballistic point contacts in 
series — without opening an excitation gap at the Fermi 
level. 

The geometry, sketched in Fig.[TJ can be seen as a gate- 



controlled realization of the puddles of metallic conduc- 
tion that may occur naturally near the QSH edge [2D - 
[25] , An electron entering the metallic area under the 
gate from one side can be either transmitted to the other 
side or reflected back to the same side, with amplitudes 
contained in the 2x2 unitary scattering matrix S(e), de- 
pendent on the energy e relative to the Fermi level. Time 
reversal symmetry requires an antisymmetric scattering 
matrix [23], S nm — —S rnn , so the reflection amplitudes 
on the diagonal are necessarily zero and the gate has no 
effect on the conductance. 

A perpendicular magnetic field B effectively removes 
this constraint, once the flux through the gate is of the 
order of a flux quantum h/e. The electronic scattering 
matrix then has the four-parameter form 



S 



r 1 t' 

t r 



^KplOa p l<P2<yz pll^y p l 03CT z 



(1) 



1 £ [0,tt/2), 0„g[O,27t), n= 1,2,3. 



We have introduced Pauli spin matrices a x , o~ z , with <t 
the 2x2 unit matrix. 

If the scattering in the quantum dot is chaotic, the 
matrix S is uniformly distributed among all 2 x 2 unitary 
matrices. The Haar measure on the unitary group gives 
the probability distribution 



P( 7 ,^ 1 ,0 2 > 3 ) = (27r)- 3 sin2 7 , 



(2) 



representing the circular unitary ensemble (CUE) of 
random-matrix theory |25j . This produces a transmis- 
sion probability T — \t\ 2 = sin 7 that is uniformly dis- 
tributed between zero and one [26j [27]. Different real- 
izations of the ensemble, with different T £ [0, 1], can be 
reached by varying the gate voltage, so that the quantum 
dot in a magnetic field functions as a tunable transmitter 
for the helical edge channels. 

We now use this quantum dot as an energy-sensitive 
detector of the presence of a Majorana zero-mode at the 
interface with a superconductor. To explain how the en- 
ergy sensitivity appears, we follow the usual procedure 
[25] of combining the electronic scattering matrix 5(e), 
the hole scattering matrix S*(— e), and the Andreev re- 
flection matrix 



''A 



'yi 



= V / l-(e/Ao) 2 + z£/A c 



(3) 



The Pauli matrix t v acts on the electron-hole degree of 
freedom and Ao is the superconducting gap. An electron 
incident on the quantum dot along a helical edge state is 
reflected back as a hole with probability 



Rhe{s) 



T(e)T(-e) 



II- 



a 2 (e)r(e)r*(~e)\< 



(4) 



At the Fermi level e = one has a — 1 and rr* = 1 — T, 
hence Rhe = 1 irrespective of the transmission probabil- 
ity T through the quantum dot. This is the Majorana 
resonance [28 a . Away from the Fermi level the resonance 



has (for T<l)a Lorentzian decay ex [1 + (e/T) 2 ]^ 1 , of 
width T = T<5dot/47r set by the average level spacing 8& t 
of the quantum dot. 

The differential conductance G = dl/dV, at bias volt- 
ages \V\ < Ao/e and in the zero-temperature limit, di- 
rectly measures the probability Q: 



G/G a = 2 + 2R he (eV), G = e 2 /h. 



(5) 



The two contributions to the conductance correspond to 
the two edges connecting the normal and superconduct- 
ing contact: The edge containing the quantum dot con- 
tributes 2e 2 /hxRhe, while the other edge remains unper- 
turbed and contributes the full 2e 2 /h — for sufficiently 
small B that the helical edge state remains gapless. 

The ensemble averaged conductance (G) has a peak 
value of 4Go at V = 0, above an off-resonant base- 
line Gbaso that we calculate as follows. We may assume 
(5dot <C Ao, so we keep a = 1. We treat the off-resonant 
scattering amplitudes at ±e as statistically independent 
random variables in the CUE, distributed according to 
Eq. pj. Substitution of the parameterization (HI into 
Eq. ffi gives, upon averaging, 

,7T/2 ,tt/2 

Gbaso/Go = 2 + 21 e?7 + / cfy_ sin27 + sin27_ 
Jo Jo 

2jt ji -2 -2 

d(f> sin 7_|_ sin 7_ 



27r |1 — cos7 + cos7_e^| 2 



§tt 2 -4«2.58. 



(6) 



A similar calculation gives the triangular line shape of 
(G(V)) as an average over the Lorentzian line shape of 
G(V), 



(G(V)) - G base ex / dT [1 + (4ireV/T5 dot ) 



2l-l 



= 1 - (AneV/Sdot) arctan ((5 do t/47reU) 
= 1 -2n 2 e\V\/S dot + 0(V 2 ). (7) 

To test these analytical predictions, we have performed 
numerical simulations of a model Hamiltonian for an 
InAs/GaSb quantum well \M HH HI]. Results are 
shown in Fig.[2]and fully confirm our expectations: With- 
out the quantum dot the Majorana resonance remains 
hidden in the background conductance (dashed curve in 
Fig. [2^,), demonstrating that the 0.1 T applied field is 
weak enough to cause no appreciable backscattering of 
the helical edge states. We then create a 200 nm x 200 nm 
quantum dot, as in Fig. IT] by applying a gate voltage. 
This suppresses the background conductance, revealing 
the Majorana resonance at V = (solid curves). Disor- 
der averaging removes all resonances from Andreev lev- 
els at V 7^ 0, so that the Majorana resonance stands 
out above the baseline conductance Gbasc, in very good 
agreement with the calculated value ffih. The triangular 
line shape of the average conductance is also confirmed 
by the simulations (Fig. [2]t 
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FIG. 2: a) Zero-temperature differential conductance G = 
dl /dV as a function of bias voltage V , calculated numerically 
for the system of Fig. [I] The bottom of the conduction band 
in the gated region (200 nm x 200 nm) is E c = — 1.5meV be- 
low the Fermi level. The black curve is for a single disorder 
realization, the red curve is the disorder average. The calcu- 
lated background conductance Gbase from Eq. B is indicated 
(arrow). For comparison, the conductance without the gate 
electrode is also shown (green dashed curve). The Majorana 
resonance is then fully absorbed in the background and invis- 
ible, b) Enlargement of the Majorana resonance from the left 
panel, to show the difference in line shape. 
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FIG. 3: Differential conductance, averaged over gate voltages 
(— 4.5meV < E c < — 1.5meV) for a single disorder realiza- 
tion. All system parameters are the same as in Fig. [2] except 
for the size of the gated region, which is 100 nm x 100 nm. 
The solid curve is at zero temperature and the dashed curve 
at a temperature of 0.1 K. 



The ensemble average in Fig.plis an average over disor- 
der realizations. As is well known from quantum dot ex- 
periments 32, 33 , statistically equivalent ensembles may 
be generated for a fixed disorder potential by varying the 
gate voltage, which is more practical from an experimen- 
tal point of view. In Fig. [3] we show a computer simula- 
tion performed in this way. To reduce the sensitivity to 
thermal averaging, we took a smaller (lOOnm x 100 nm) 
quantum dot, keeping the magnetic field at 0.1 T. The 
simulation shows that the Majorana resonance remains 
clearly visible above the background conductance at tem- 
peratures of 100 mK. 

So much for the detection of Majorana zero-modes. In 
the final part of this paper, we take a longer term perspec- 
tive and present a geometry that allows for the braiding 
of pairs of Majorana fermions, for the demonstration of 
the predicted non-Abelian statistics [3] . While the quan- 
tum spin Hall edge seems ideally suited for the detection 



of Majorana zero-modes, its one-dimensionality prevents 
the exchange of adjacent Majoranas. What is needed 
is a Y- or T-junction of superconductors to perform the 
"three-point turn" introduced by Alicea et al. [M] and 
implemented in a variety of braiding proposals for a net- 
work of nanowires [35H38] . 

In Fig. [4^, we show how a constriction in the quan- 
tum spin Hall insulator can be used to achieve the same 
functionality as a crossing of nanowires. The constric- 
tion couples the helical edge states on opposite edges by 
tunneling, which is effective if it is narrower than the de- 
cay length of the edge states (100 nm or smaller). The 
coupling may be increased, if needed, by gating the con- 
striction region into the conduction band. Three of the 
four edges leading into the constriction are gapped by 
a superconducting island. The fourth edge contains one 
of the quantum dots discussed earlier, tuned to a gate 
voltage interval of minimal transmission T<1. 

Let us check that the constriction traps a Majorana 
zero-mode. Helical edge states incident on the constric- 
tion from the three superconductors have reflection am- 
plitudes that are contained in a 3 x 3 reflection matrix 
r(e). Neglecting transmission throught the quantum dot, 
this is a unitary matrix. (A nonzero T will give a finite 
width to the zero mode.) A bound state in the con- 
striction at energy e is a solution of the determinant al 
equation [3!5] 



Det [1 - a 2 Ar(e)A*r*(-e)} = 0, 



(8) 



a diagonal matrix con- 



with A = diag (e^ 1 , e^ 2 , e^ 3 
taining the phase <p n of the order parameter on the n- 
th superconductor. Since a(0) = 1, see Eq. Efy, the 
condition for a zero mode is that the matrix uu*, with 
u = Ar(0), has an eigenvalue equal to +1. The eigen- 
values of uu* come in complex conjugate pairs e ±J ^. An 
unpaired eigenvalue at —1 is forbidden by Detuu* = 1, 
but an unpaired eigenvalue at +1 is allowed and in fact 
necessary when the dimensionality of u is odd — as it is 
here. 

In Fig. [4b, c we combine two constrictions in 7r-shaped 
circuit, to perform the braiding protocol of Ref. |38j . 
There are six Majorana fermions, one at each constric- 
tion and four more trapped by quantum dots along the 
quantum spin Hall edge. Adjacent Majorana operators 
lm , 7m' j for example jb and je , are coupled by the 
charging energy Eq of the intermediate superconducting 
island S n through the Hamiltonian |36j 



H., 



iUnlmlm' , U n oc exp [- \/8Ej ($ n )/Ec 



(9) 



This coupling can be switched on and off by adjusting 
the Josephson energy Ej = E cos(e$/h) of the super- 
conducting island, via the magnetic flux $ through a 
split Josephson junction that connects the island to a 
superconducting ground. 

As worked out in Refs. [35l [36], the alternating cou- 
pling and decoupling of adjacent Majoranas has the ef- 
fect of exchanging them: One effectively moves the Ma- 
jorana at B through the T-junction towards F, followed 




V 








^1 





gate I 



J gate 



• o 

braid (twice) 



• ■--<? — • • 




Z2 



...measure 



6 



FIG. 4: a^ Constriction in a quantum spin-Hall insulator, contacted along three edges by a superconducting island. If the 
fourth edge is blocked by a gate electrode, the constriction traps a Majorana fermion. Each superconducting island is connected 
to a superconducting ground by a split Josephson junction enclosing a magnetic flux, indicated schematically by a x symbol 
and shown expanded for one of the islands, b) Two constrictions in series form a 7t-shaped circuit that can be used to braid the 
Majorana fermions (green circles). The flux through each split Josephson junction controls the coupling of adjacent Majoranas. 
c) Schematic of the braiding operation in the 7r-circuit [33]. Coupled Majoranas are connected by a solid line, uncoupled 
Majoranas by a dashed line. 



by C i— > B and finally F >->• C completes the exchange 
of B and C. If this exchange is repeated, one ends up 
with the original configuration of Majoranas, but in an 
orthogonal state: The fermion parity of S4 has switched 
between even and odd. This signature of non-Abelian 
statistics can be measured as described in Ref. [38], as 
a shift in the resonance frequency of a superconducting 
transmission line containing the circuit. 

In conclusion, we have shown how the helical edge 
state in a quantum spin Hall insulator may be used as 
a single-channel, disorder-insensitive alternative to semi- 
conductor nanowires, for the detection and braiding of 
Majorana fermions. For all we know, the experiments 
on InAs/GaSb quantum wells [10] may already have pro- 
duced the predicted Majorana zero- modes [T], but since 
the 4e 2 /h conductance resonance is hidden in the 4e 2 //i 
off-resonant background there is no way to tell. The 
quantum dot geometry proposed here lowers the aver- 



age background to about 2.6e 2 //i, allowing for the emer- 
gence of the Majorana resonance. This seems to be an 
experiment that is fully within reach of existing devices, 
requiring only the addition of a nanostructured gate elec- 
trode and the application of a weak magnetic field. As a 
longer-term perspective, we have shown how a constric- 
tion in the quantum spin Hall insulator can reproduce 
the functionality of a nanowire T-junction, required for 
braiding and for the demonstration of non-Abelian statis- 
tics. 

The numerical calculations were performed using the 
kwant package developed by A. R. Akhmerov, C. W. 
Groth, X. Waintal, and M. Wimmer. We thank C. 
Liu for providing us with the model parameters of the 
InAs/GaSb quantum well. This work was supported by 
the Dutch Science Foundation NWO/FOM, by an ERC 
Advanced Investigator Grant, and by the China Scholar- 
ship Council. 
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Appendix A: Description of the numerical 
simulations 



Our numerical simulations are based on the four-band 
Hamiltonian of an InAs/GaSb quantum well [HI H2] . 
which in zero magnetic field takes the form 



H(k) 



H (k) 
-Hf(-k) 



Ht(k) 
HS(-k) 



(Al) 



as a function of wave vector k — (k x , k y ) in the x-y plane 
of the semiconductor layers. The block structure refers 
to the spin degree of freedom, while each block itself has 



&o = 


- 0.0078 eV 


CO = 


0.0002 eV 


6i = 


- 5.8 eV-A 2 


Cl = 


0.00066 eV-A 


b 2 = 


-66.0 eV-A 2 


C2 = 


0.0006 eV-A 


bz = 


0.37 eV-A 


C-A = 


-0.07 eV-A 



TABLE I: Parameters of the four-band Hamiltonian (Al) , 
representative for a heterostructure consisting of 10 nm InAs 
and 10 nm GaSb layers, sandwiched between AlSb barriers 



a 2 x 2 matrix structure that refers to the (s,p) orbital 
degree of freedom. 

The diagonal block describes the hybridization of the 
s and p orbitals, 



Ho 



U + b Q - b + k 2 b 3 k + 

b 3 k_ U-b -b-k 2 



(A2) 



fcn 



k x ± ik v , and b± = b\ ± 62 ■ We 



with k 2 — k 2 + k 2 , k± = k x zc iK y , ana o± — ui zn 02- 
have included an electrostatic potential U, to account 
for the effects of a gate electrode. The off-diagonal block 
describes the spin-orbit coupling by inversion asymmetry 
(Rashba and Dresselhaus effects), 



Si(fc) 



c\k + + ic 3 k_ —Co 
c c 2 fc_ 



(A3) 



to first order in k. The parameter values we used in our 
simulations, taken from Ref. UHl are listed in Table I. 
Time reversal symmetry is expressed by 



H(k) = cj y n* (-k)a 



v 



(A4) 



where the Pauli matrix a y acts on the spin blocks. A 
perpendicular magnetic field B — (0, 0, B) breaks time 
reversal symmetry, via an orbital and a Zeeman effect. 
The orbital effect is accounted for by the substitution 
k n- k — (e/h)A, with vector potential A = (0,Bx,0). 
The Zeeman energy g^B<r z is a negligibly small effect 
in the weak magnetic fields B ss 0.1 T considered here, 
so we do not include it. 

The effect of a superconducting contact is introduced 
via the Bogoliubov-De Gennes Hamiltonian, 



H 



BdG 



H[k 



- (e/H)A\ 

At 



-a y n*[-k- {e/h)A}cjy 



(A5) 
The blocks of if BdG refer to the (e, h) electron- hole degree 
of freedom, coupled by the pair potential A induced by 
the superconducting contact. Electron- hole symmetry is 
expressed by 

H B dG{k) = -(<Jy <8> T y )H BdG (-k)(a y ® t v ), (A6) 

where the Pauli matrix r y acts on the electron-hole 
blocks. A spin-singlet s-wave proximity effect may still 
couple the s and p orbitals of the quantum well |30j , with- 
out breaking either electron-hole or time-reversal symme- 
try, but for simplicity here we take a scalar A. 



For the numerical simulations we discretize the Hamil- 
tonian ( |A5[ ) on a square lattice (lattice constant a — 
10 nm) in the geometry shown to scale in Fig. [Tl We set 
A = 1 mcV in the superconducting contact and zero else- 
where. The effect of the gate electrode is modeled by an 
offset [/gate of the electrostatic potential in the area under 
the gate. Disorder in the quantum well is modeled on the 
lattice by a random on-site potential [/disorder, uniformly 
distributed in the interval (—5 meV, +5 meV). The mag- 
netic field in the quantum well and the normal-metal con- 
tact is fixed at B = 0.103 T, corresponding to one flux 
quantum h/e through an area of size 200 nm x 200 nm. 
Inside the superconducting contact we set B = 0, ignor- 
ing the penetration of flux in magnetic vortices. 

At excitation energies |e| < A the electrons and holes 
incident on the superconductor are fully reflected. The 



reflection amplitudes are contained in an 8 x 8 unitary re- 
flection matrix. The 4x4 Andreev reflection block The (e) 
gives the differential conductance at zero temperature, 

G(V,0)= 2 ^-Trr he (eV)r{ e (eV). (A7) 

The corresponding result at finite temperature T Q follows 
upon integration, 



G[V.T„) = - j deG(e/e,0)^f, 

-oo *-^ 



(A8) 



/(e,V,T ) 



l + exp[(e-eV)/AfeT ]' 



